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Abstract 

The influence of a fullerene molecule trapped inside a single-wall carbon nanotube on resonant 
electron transport at low temperatures and strong polaronic coupling is theoretically discussed. 
Strong peak to peak fluctuations and anomalous temperature behavior of conductance amplitudes 
are predicted and investigated. The influence of the chiral properties of carbon nanotubes on 
transport is also studied. 

PACS numbers: 73.63.Fg, 73.23.-b, 73.63.-b 



I. INTRODUCTION 



Carbon nano-peapods were discovered almost one decade ago [l] (see also the review 



ano-peapc 

B,S,B). 



papers Refs. [2j, [3|, |4J]). They have been proven to be novel graphitic structures with elec- 
trical and mechanical properties that can not be reduced to the sum of properties of their 
subcomponents (an empty carbon nanotube and individual fullerene molecules). Peapods, 
a composite system consisting of a rigid single- wall nanotube (SWNT) cage and an inner 
mechanically soft chain of fullerenes, can exhibit interesting nano-electro-mechanical prop- 
erties if the coupling of delocalized electron states on the tube with the fullerene molecular 
orbitals is sufficiently strong. Buckyballs inside a SWNT are neutral objects and one can 

not expect any direct electrical influence of encapsulated molecules on electron transport 

n 

along a nanotube shell. However, STM spectroscopy of semiconducting peapods |5( showed 
that the local electron density of states in carbon nanotubes is strongly modified by encap- 
sulated Cqo molecules. These measurements can be theoretically interpreted in terms of a 
strong hybridization of electronic states on the nanotube surface with the lowest unoccupied 
molecnlar orbitals (LUMO) of the encaged C m . The corresponding hybridization strength 
was estimated in Refs. [5|, [6[ to be close to 1 eV. This energy is comparable to other charac- 
teristic energies of SWNTs and this means that virtual hopping of delocalized electrons from 
the SWNT to the local LUMO states is possible and can strongly influence the electronic 
properties of the SWNTs. Even though all samples in the cited experiment were semicon- 
ducting nanotubes, one could expect from theoretical considerations an analogous strong 
hybridization at least for chiral metallic nanotubes (see Ref. 6j where a theory of scanning 
tunneling spectroscopy of long C 60 peapods was developed). Notice that such a strong hy- 
bridization between fullerene-derived and nanotube-derived energy levels in peapods could 
shift the LUMO states even below the Fermi energy of the SWNT. This possibility was 
however ruled out in recent photoemission experiments, {7]. 

In this work we suggest that a strong coupling between the electronic states on the 
SWNT shell and the localized states on the fullerenes affects resonant tunneling transport of 
electrons through the metallic peapod system in a manner that is experimentally observable. 
In particular, we propose that resonant electron transport due to polaronic effects provides 
a method for studying the fullerene nanomechanical dynamics through the temperature 
dependence of the resonant conductance peaks. This should be readily observable since the 
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characteristic temperature scale for the strong polaronic modification of the peaks is set 
by the fullerene vibration frequencies and may be very low compared to typical electronic 
energies. 

Recent experiments with high-quality peapods demonstrated Coulomb blockade and 
Kondo physics at low temperatures, j^l, 9]. The spacings of conductance peaks was found 
to be regular over a wide range of gate voltages, indicating that the Constant Interaction 
Model of Coulomb blockade (see e.g. Ref. [l0|) applies for these dots. The peak amplitudes, 
on the other hand, strongly fluctuate from peak to peak. In the case of an empty metallic 
SWNT both peak spacings and amplitudes should ideally vary slowly with the gate voltage, 
except for variations of level spacings dictated by the particular SWNT bandstructure which 
leads to four-electron shells. While earlier experiments on empty nanotubes have exhibited 
strongly fluctuating peak conductances (see e.g., Ref. 11 
shown remarkably gate-independent peak amplitudes [12 



, defect-free SWNTs have indeed 
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We show here that for carbon nano-peapods, electromechanical effects can lead to strong 
fluctuations of the Coulomb peaks. The mechanism will be effective even with defect-free 
nanotube shells. This phenomenon may explain the strongly fluctuating peak amplitudes 
seen in peapod experiments 8|, |9] - detailed studies of Coulomb blockade peak fluctuations 
could thus reveal polaronic effects in quantum dots. Molecular conductors constitute another 
important class for which the vibronic coupling is also nearly always present, and for which it 
plays a fundamental role in the study of the transport properties, see 14], [15] and references 
therein. 

The paper is organised as follows. In the following Section [Til we investigate how the 
presence of even one encapsulated fullerene molecule, treated as a short-range scattering 
potential, can influence the electron energy spectrum of the SWNT by a strong polaronic 
effect. This will allow us to define the Hamiltonian of the peapod in an independent electron 
model. In Section IHIl the current through the peapod and hence its conductance are evalu- 
ated and studied, in this regime of strong polaronic coupling. Since the fullerene scattering 
potential in reality has an intermediate range, we discuss in Section IIVI how the peapod 
conductance is affected in the opposite limit of a long-range fullerene-induced scattering po- 
tential. Here, the appropriate model for the SWNT electrons involves the use of the Dirac 
equation. Finally, in Section V we present our conclusions. 
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II. ENERGY SPECTRUM IN PRESENCE OF A SHORT-RANGE SCATTERING 
POTENTIAL: THE MODEL 



The position of an individual buckyball inside the tube, for a partially filled peapod, 
is not fixed, and mechanical motion of the fullerene is excited in the process of virtual 
electron transitions onto and off the fullerene LUMO level. The effective excitation energy 
huo strongly depends on the buckyball confinement potential. For vibrations in the nearly 
harmonic (for small displacements) transverse Girifalco potential [if], 1^] the corresponding 
energy fkvt is of the order of meV's. Notice that electro-mechanical vibrations of Cqq in 
the van der Waals-like potential between two gold electrodes have already been measured 
(fko ~ 5 meV) in tunneling experiments is]]. The longitudinal motion of Cqq inside the 
tube is less restricted and the corresponding excitation energy ftw\ could be much smaller. 
It is worth to notice here that the nanotube shell of peapods can in fact be quite full of 
large (~ 1 nm) holes in the tube walls. Only a small number of openings, used before 
to fill the tube, are healed in the post processing, [3]. The average spacing between the 
holes can be estimated to be ~ 5 — 10 nm, The holes produce a repulsive potential for 
fullerenes inside the tube, which we will model in our analytical calculations as a harmonic 
confinement potential with a characteristic energy fru)o ~ 0.1 meV, corresponding roughly 
to the zero-point energy of a Ceo confined to a length equal to the estimated average spacing 
between holes. 

We start with a single encapsulated Ceo placed at a distance I from the nanotube end. 
The virtual processes of electron hopping onto and off the fullerene LUMO level (shifted 
upwards by the charging energy ~ e 2 /2d, where d is the fullerene diameter) produce a local 
attractive potential for conduction electrons on the tube, The first question one has to 
consider in order to formulate a theoretical model for the problem is whether this potential 
is short- or long-range with respect to its effect on electrons in metallic nanotubes. If the 
characteristic length of the scattering potential due to the fullerene is defined by its diameter, 
of the order of 1 nm, its range is determined by a comparison with the Fermi wavelength, 
Ap, of electrons in a metallic SWNT. Since Ap ~ 7 A, [20j, and hence is remarkable close to 
the fullerene diameter, the range of the potential is of intermediate range. It is therefore not 
a priori clear which of the two different limits, where an analytical investigation is possible, 
is the more appropriate one and we shall therefore discuss both. 
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In the first limit, the scattering potential is described as a short-range potential, which 
means that the hybridization-induced scattering potential is able to strongly backscatter 
electrons with a large momentum transfer, 5k ~ 2k F . This is the regime studied in this 
Section by means of a single-electron model. In the opposite limit the scattering potential 
is described as a long-range potential when the real band structure of the SWNT becomes 
important. This limit is hence discussed in Section HVl using the Dirac equation. 

A short-range hybridization-induced scattering potential with the hybridization strength 
th — 1 eV changes the distribution of level spacings in the conduction band of a finite-sized 
nanotube. By solving the Schrodinger equation on a finite interval (L) with a 5-function 
scatterer placed at point x = I one gets the spectrum equation 

k sin kL = U{k) sin kl sin k(L — I), (1) 

where k is the electron wave vector, (e = ti 2 k 2 /2m), and the amplitude of the local scattering 



potential in our case represents the hybridization potential of Ref. 



6|. In our notation 



U(k) = ( k f) V h (e) , V h (e) = -M*- . (2) 

Here, d ~ Af and Eq > e is the energy of the hybridized side-level (Eq ~ 2s-f if one takes 
into account the charging energy for electron hopping to the Cgo LUMO level). 

The energy spectrum of Eqs. ([T]) and (j2j) when l/L is not a small rational number is 
irregular for energies not far from the Fermi energy (see FIG. 1), and indicates certain 
deviations from the mean level spacing Ax, = Trhv-p/L (notice that due to the two identical 
conduction bands in a SWNT its mean level spacing is Ax/2). However, as one can see from 
FIG. 1, for a single scatterer the "randomization" is not pronounced even for energies close 
to the Fermi energy. About 40% of the levels deviate from the mean level spacing less than 
18%, and among these levels 75% deviate less than 12%. 

A shift 5x in the position of a scatterer leads to shifts of all quantized energy levels 

dk 

5e n = hv F —^Sx, (3) 
ox 

where the derivative of electron momentum with respect to scatterer position can easily be 
found from Eq. ([1]). 

In order to get analytical results we will model the potential that confines the longitudinal 
motion of our scatterer (buckyball inside a tube) as a harmonic potential. In the adiabatic 
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FIG. 1: Energy levels, e n L/hvF, for a peapod with just one fullerene inside. Unlike an empty 
SWNT, the spectrum is not uniform; still, for a single scatterer (fullerene), the randomization is 
not pronounced. The Fermi energy Ep is such that EpL/hv-p ~ 90. 

limit, Uq <^vp/L, when the electron spectrum is well denned for any position of the vibrating 
impurity, the total Hamiltonian for a two-terminal conduction geometry reads 1 : 

H = H L + H R + H T ^ L + H TR + Hqd (4) 

where H a = J2k a £ k a c \ a c k a is the Hamiltonian for the leads and a = L, R for left and right 
reservoir respectively. The quantum dot (QD) Hamiltonian reads: 

Hqb = E ^n°n + Yl + b ) + ^\ (5) 

n n 

where e n is the set of levels shown in FIG. 1, hw is the vibration energy quantum, V n (L, I) = 
HvFXo(dk n /dl) and xq is the amplitude of zero-point fluctuations of the bosonic field (xq ~ 
0.2 A for C*6o and for an energy value ftwo ~ 0.1 meV); c n (cjj and b(tf) are fermionic and 
bosonic operators with canonical commutation relations. The tunneling Hamiltonian finally 
reads: 

H ^ = Yj(c\ a c n + Yi. C .) , (6) 

where we suppose that the hopping matrix elements t between the two leads and the dot 
are of the same order. 

[1] The total Hamiltonian H in Eq. [4] does not take into account the scattering potential represented by the 
holes on the tube shell. The possible consequences of such a issue are discussed in the conclusions 
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We should point out that the coupling to the electronic states of the SWNT does not 
significantly affect fullerene vibrations, which are mostly determined by fullerene confinement 
by structural defects on the carbon nanotube. This is true even though the coupling to each 
single electronic state may be strong and depend on the fullerene position. However, when 
evaluating the total fullerene energy shift caused by coupling to a very large number of 
different electronic states, the coordinate dependence averages out and such a coupling does 
not contribute to the mechanical force. 

The Hamiltonian -£/qd can be diagonalized by means of a unitary transformation, see 
e.g. [21]: H QD = e s H QB e~ s , with S = /)/^o)4c n (6 t - b). The result of this 

transformation is 

#QD = J2( £n ~ A ) C « C " + kuJ ^ h > ( 7 ) 
n 

where A = X^Fiujq is the polaronic shift of the resonant energy level and A n = V n (L, l)/(huj ) 
is the electron- vibron coupling constant. The tunneling Hamiltonian is also transformed by 
the unitary trasnformation and it reads 

#T,a = X>( C Lcn* + H.C.) (8) 

where the operator X is defined as 

V n (L,l) 



X = exp 



III. RESONANT TUNNELING 



(6 f - b) 



(9) 



Resonant electron tunneling through a vibrating quantum dot was considered in many 



papers (see e.g. Refs. 22|, |23| and references therein). For strong electron-phonon interac- 
tion, two non-perturbative effects determine the low-temperature electron transport - (i) a 
polaronic shift of the resonant energy level, e n — > e n — \ 2 huJo and (ii) a "polaronic blockade" 
(exponential suppression) of the peak conductance at low temperatures, G\ cx exp (—A 2 ). 
In our case, the coupling "constant" A n is a level-dependent quantity and it strongly fluc- 
tuates from level to level. Such fluctuations are well visible in FIG. 2. We observe that 
these coupling constants being proportional to the derivative of the energy spectrum exhibit 
much more pronounced fluctuations then the energies themselves. The mean value of A 2 
for strong backscattering potential (U(kp) ~ Ap 1 ) can be estimated as A 2 ~ {syXq / hoJoL) 2 . 
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FIG. 2: Behaviour of the polaronic coupling constant A^ as a function of state index n. The coupling 
constant strongly fluctuates since it is proportional to the derivative of the energy spectrum. The 
corresponding value for the 8-th level, of the order of unity, is off the scale. The Fermi level is 
between the 10-th and the 11-th levels. We stress that we do not intend to show the actual and 
absolute numerical values of the couplings A. These "constants" can be obtained by fitting the 



experimental data to theory, as in [2j|. FIG. 2 only aims to illustrate the strong level-to-level 
fluctuations of the coupling constant. In the case of N independent molecules in the peapod, the 
average value of the coupling constant can be iV times larger than the values shown in the plot. 
In particular, the values of the coupling constant A shown in the plot have been evaluated for a 
sample whose length is ten nanometres, and with a LUMO energy and a hybridization energy of 
the same order as the Fermi energy. These are the same values used to evaluate the spectrum in 
FIG. 1. 

So, for the case of N independent scatterers N = L/d c ,( d c is the characteristic range of 
the longitudinal confinement potential) the total polaronic shift will scale as 1/L. Since 
the charging energy of our ID system also scales as 1/L, it will be difficult to single out 
polaronic effects in the distribution of peak spacings as a function of gate voltage. Besides, 
the irregular character of level spacings could be caused by different physical effects, e.g., 
by electron scattering due to imperfections (holes) of the nanotube surface. Is there another 
manifestation of low-energy electro-mechanical effects in resonant electron tunneling? 
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A. Electric current 



In order to study the peak to peak fluctuations in the electrical conductance, we first 
recover in this section the electric current flowing through the system. 
We start from the well known Meir-Wingreen formula, 25j, |26|, 



IiV) = -^jde [ft(s) - f R (s)]Qm{Tr[TG^(e)]} 



(10) 



where /l/r(£) — [exp(e — e F =)= eV/2) + l] -1 are the Fermi functions for the left and right 
reservoirs, Gq^e) is the retarded Green function of the dot and T = r L r R /(r L + T R ). 
Tl/r = 27rt 2 pL/R are the widths of the resonance due to tunneling from the left or right lead, 
and pl/r is the density of the states (DoS) in the reservoirs. 

In order to evaluate the current in Eq. (fTOl) . one needs the Green function of the dot 
which takes into account, at least in principle, not only all the possible processes in the dot, 
as inelastic processes anc multiple scattering processes, but also the effects of the vibrons 

251 ] . To get analytical results, we will treat the leads, as in [23J, as 



on the reservoirs, 



23 



unaffected by the vibrons modes in the dot. Besides, we suppose that the bandwidth in the 
contacts is much larger than both the resonance width and the vibron energy quantum, that 
is we will evaluate the electron part of Green function in the so called wide-band limit, j^. 

Under such assumptions, the retarded single particle Green function can easily be evalu- 
ated, see also [2l|], and it reads 

.r L + r R 



G&(f) ~ -*0(t)exp 



A 



x 



exp [-\ 2 (2n(f3hLO ) + 1 

t a 



h 



(11) 



J2 Ii{2X 2 [n(/3fku )(n(/3tku ) + l)] 1/2 }exp 

£=— oo 

where n(x) = (e x — 1) _1 is the Bose distribution, (3 = l/k^T, l\ is the modified Bessel 
function of the first kind, and for sake of simplicity, we have omitted the index n from the 
coupling constant A. Equation (TTT1) can be easily Fourier transformed and the expression of 
the current written as in 23 ] 



-A 2 (2n(/3^,,)+l) 



I(V) = de[f L (e) - Ms)]e 

oo 

x ^{2A 2 [n(/5^o)(n(/?^o) + l)] 1/2 }e-^ o/3/2 T B w(5,/) 



(12) 



l=— oo 
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where 



with Ap i = A — e n . 



r L r R 



T B w{e, I) = — o , (13) 

(e + ^ + A pol )2 + (^) 2 ' 



B. Electric conductance 



A simple expression of the conductance can be derived from Eq. ffl2|) in the linear 
regime, eV — > 0, and weak coupling to the leads, IY + Tr <C min{T, frwo}. In addition 
to the latter condition, in order to have polaronic states, it is necessary, physically, that 
r < T < X 2 hu>Q, as well. These conditions state that the time spent by the electron in 
the dot must be larger than the characteristic time needed to form a polaron. Under these 
assumptions, the Lorentzian in Eq. fflBl is just a 5-function centered at e = e n — A — ItiujQ. 
Then, the integral in Eq. (TT2l) can be easily evaluated and the current reads 

I{V) = 27re^[/ L (£„ - A - IHuj q ) - f R (s n - A - lfuv )) 

oo 

x e - A2(2n(/3 ^ o)+1) Ii{2\ 2 [n((3huj ){n{[3huJo) + i)] 1 /*^-"^/* . ( 14 ) 

l=— oo 

Since we are interested in the conductance peak height, we consider the physical condition 
where, by means of the gate voltage, the levels in the dot are tuned such that Ef = e n — A. 
Then, in the limit eV — > 0, a simple analytical expression can be written: 

G X (T) = G x=0 (T)F x (hw /k B T), (15) 

where G\=o{T) ~ (jr /2)GqT jk^T is the standard resonance conductance of a single-level 
quantum dot (Go = e 2 /h is the conductance quantum) and the function F\(x), with x = 
hu> /k B T, is represented as a series 

rp ( \ r a2m , 9 nil V exp(-/x/2)/,[2AVn(x)(l + n(x))] 

F x (x) = exp{-A [1 + 2n(x)]} ^ . (16) 

l=— oo v ' ' 

It is easy to check that F A=0 (x) = 1. 

For strong electron- vibron coupling (A > 1), polaronic effects can significantly suppress 
the conductance at low temperatures, (r k-gT TujJo), for which one finds F(x) ~ 
exp(— A 2 ). In this regime, the usual ~ 1/T behaviour is recovered for the conductance, but 
the latter is strongly suppressed because of the factor exp(— A 2 ). This is the signal of the 
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FIG. 3: (Color online) Behavior of the conductance G\(T) = G\tkvo/T as a function of the 
renormalized temperature T = k-QT/hujQ and for three different values of the coupling constant A. 
An anomalous temperature dependence of the resonance peak height is well visible in the region 
k-^T < X 2 fkOQ. Inset: the behaviour of the conductance for strong coupling, A = 3, is shown in 
a narrow region close to zero temperature. The behaviour is dominated by exp(— A 2 )/T close to 
zero. Then, the dominant term is proportional to ~ (1/y/T) exp(— \ 2 hujQ/ Ak^T). 

reduced probability for an electron to tunnel from a bare electronic state in the reservoir 
to a polaronic state in the dot. As soon as the temperature becomes of the order of Hluq, 
thermally excited vibronic modes of the fullerene molecule appear. The characteristic energy 
of fullerene thermal vibrations is n(f3huQ)huQ, n{x) being the Bose distribution. This energy 
should be compared with polaronic energy shift \ 2 hu)o induced by fullerene displacement. 
The corresponding ratio determines the relative strength of destructive influence of thermal 
vibrations on polaronic effects. Indeed, at temperatures hw <C k B T < X 2 Jtwq, by keeping 
the first 1- independent term in the asymptotic expansion of the Bessel function J;, F(x) ~ 
((yV/x )/A) exp(— A 2 a;/4), and the conductance reads, 




between the vibrons and the polarons, giving rise to a non-monotonic behaviour. At tern- 





(17) 



The latter equation indicates that for temperatures larger than Hljq, a competition arises 
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peratures k^T ^> \ 2 huo, when the polaronic blockade is completely lifted, the conductance 
scales as 1/T, as can be numerically verified. At such high temperatures, the polaronic 
effects get completely destroyed by the thermally excited vibrons. The temperature lift- 
ing of the polaronic blockade results in a non-monotonic temperature dependence of the 
conductance peak. This behavior is shown in FIG. 3, where we plot the temperature de- 
pendence of the amplitude of the Coulomb blockade conductance peak for different values 
of the electron- vibron coupling A 2 . It is clearly seen from the figure that polaronic effects 
are pronounced only for strong electron- vibron interaction (A > 1), and in the temperature 
region kv,T < \ 2 hjo anomalous temperature dependence of resonance peaks appears. The 
observation in an experiment of a low-temperature scaling which differs from the usual In- 
dependence for narrow Coulomb blockade conductance peaks, in general, indicates that the 
level width is not the only relevant energy scale in the problem. An anomalous temperature 
behavior and l/\/T-scaling is the signature of strong polaronic effects, as experimentally 



observed, [24]. 



IV. POLARONIC EFFECTS IN LONG-RANGE SCATTERING LIMIT 

The polaronic effects discussed in the previous Sections appear due to strong backscat- 
tering of the SWNT electrons by the fullerene-induced potential, which was assumed to be 
short-range. Here we want to consider the opposite limit of a long-range potential, where 
we can also obtain an analytical solution. In this case the real band structure of the SWNT 
electrons must be considered. It is known that the specific band structure (Dirac-like spec- 
trum) leads to a new type of backscattering by even a long-range potential. In this Section 
we will study the influence on polaronic effects of this type of backscattering, which is closely 
related to the chiral properties of the carbon nanotube (the chiral angle). 

The band structure of graphene in the vicinity of the Fermi energy (these points in 
the Brillouin zone of 2D graphite are usually labelled K and K') is described by the Dirac 
Hamiltonian for massless fermions H-q = hv^ap, 28j (see Refs. 29| where a Dirac-like energy 
spectrum in graphene was confirmed in experiments). Here, p is the 2D electron momentum, 
a are the Pauli 2x2 matrices which describe the chiral properties (pseudospin) of particles 
and holes in each conduction band. The chiralities opj\p\ = ±1 of particles and holes are 
opposite in the same band and particles (holes) in different bands are characterized also 
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by opposite chiralities. Conservation of chirality in scattering by a scalar potential leads 
to various unusual phenomena in graphene that are currently discussed in the literature 
(see, e.g., Ref. [30]). The wrapping of a graphene sheet into a SWNT results in a ID 
relativistic form of the electron dispersion, which depending on the chiral indeces of the 
nanotube, describes either massless (metallic SWNT) or massive (semiconducting SWNT) 
Dirac fermions e(q) = -±.hv-F\Jq 2 + A^, where A m oc I/Rjsit is the gap in electron spectrum, 
and -Rnt is the radius of the nanotube. 

Now, we will treat the hybridization-induced scattering potential as a long-range po- 
tential p and reconsider the problem of resonant electron tunneling using the Dirac-like 
electron spectrum. Scattering of electrons by a long-range potential (5k <C kp) leaves the 
quasiparticles in the same band and thus does not mix particles (or holes) with opposite 
chiralities. In this case the chiral properties of electrons in a SWNT play a significant role 
and should be taken into account. Using the same approach as in the Ref. [6(, we consider 
the spectrum problem for a short peapod and discuss possible polaronic effects in resonant 
electron transport through this system. 

At the boundaries of a nanotube, x = 0,L, the scattering potential is sharp and this 
results in strong backscattering 5k ~ 2kp by the boundaries. The spinor wave function of 
an electron in a nanotube now takes the form ty(x) = e lkFX ^ + (x) + e _lfcpX \l/__(x) , where 
subindices (+,-) label two identical bands of the energy spectrum. The general boundary 
condition for a finite-sized system of Dirac fermions is the absence of current through the 
boundaries. This condition is satisfied in our case if we assume perfect reflection at the 
boundaries, L *-> R T , where -R^L 1 * 1 ) denote right (left)-moving particles in the "+" or "-" 
band. 

In reference Jfi], the scattering matrix for a long-range potential was calculated using the 
Dirac equation. Now, we would like to incorporate this result in order to consider resonant 
tunneling through the one dimensional wire. This can be easily done in the limit when the 
distances from the scatterer to the ends of the wire are much bigger than the spatial extension 
of the potential. Hence, if one is not interested in the exact position of the resonant level, 
then the resonant transmission through the given level can be obtained in the approximation 
where the scatterer can be treated as a point-like defect. Such a defect will be characterized 
by the same scattering matrix obtained for the electrons propagating in an infinite wire as 
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The local (^-function) potential at x = I now represents a long-range imperfection which 
scatters fermions only within its own band. The spectral problem can be solved by the 
standard method of finding the solution of the Dirac equation in the ranges x ^ I and 
matching them at x = I by using the properties of the Dirac equation in a 5-function poten- 
tial. In the general case of massive Dirac particles, the spectrum equation is cumbersome 
and lengthy. We are interested in transport through a metallic SWNT, that is in the case 
of a massless Dirac spectrum, A m = 0. In this limit, e = ±hvpq, the spectrum equation can 
be represented in the simple form 

t{9, if) cos{2k F L) + r{9, ip) cos[2g(L - 2/)] - cos(2gL - ip) = 0, (18) 

where t(9,tp),r(9,ip) are the transmission and reflection coefficients for the ^-function po- 
tential scattering problem by massless chiral particles. In this case, 

r{9,cp) = sin 2 # sin 2 | , (19) 

where 9 is the chiral angle of the nanotube and 



We see from Eqs. f ll8l) - (l20l) that for a massless Dirac spectrum the reflection coefficient 
vanishes for 9 = 0, that is for armchair nanotubes. The corresponding spectrum is uniform 
with level spacing nhvF/2L, and the total energy shift, hvpip/2L, is induced by forward 
electron scattering by the long-range potential. The absence of electron backscattering in 
non-chiral SWNTs (formally due to conservation of chirality) has repeatedly been discussed 
in the literature (see e.g. Ref. [3(J]). In general, the suppression of backscattering in metallic 
SWNTs as compared to semiconducting nanotubes (where the presence of a gap in the 
spectrum mixes states with opposite chiralities) is often used to explain the remarkably 



good conducting properties of long metallic carbon nanotubes, [311 ]. 

For the spectrum in Eq. f|T8|) . one can easily estimate the electron- vibron coupling as 
A ~ r(9,ip)vFX /u L 2 . Polaronic effects are determined by the square of this parameter, 
which is smaller by a factor (Xp/L) 2 then the analogous coupling expected for a short-range 
hybridization potential. 

The origin of such a weak polaronic coupling constant resides in the strongly reduced 
backscattering processes, since polaronic effects are extremely sensible to that. Then, in 
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transport measurements, where electronic backscattering with a large momentum transfer 
is expected because of the same order of magnitude of the electron Fermi wavelength and 
the fullerene diameter, the description of the hybridization-induced scattering potential as 
a long-range potential is not adequate. In this case, short-range approximation (sections [TT] 
and IHip is more appropriate for description of polaronic effects. 



V. CONCLUSIONS 



In conclusion, we have showed that the elastically soft subcomponent of peapods (fullerene 
molecules trapped inside a SWNT) can strongly influence low-temperature properties of res- 
onant electron transport through short metallic peapods. If encapsulated Ceo 's do modify 
the electronic structure of metallic nanotubes, (the effect was already observed for semicon- 
ducting peapods [51]), the mechanical degree of freedom associated with fullerene dynamics 
will influence, via polaronic effects, both the distribution of conductance peak amplitudes 
and peak spacings as a function of gate voltage. In our model we considered independent vi- 
brations of individual fullerenes (a disordered chain of C 60 's inside a nanotube). This model 
predicts strong fluctuations and anomalous low-temperature behavior of Coulomb blockade 
conductance peaks. Such fluctuations have been experimentally observed, 24 1. 

One more question remains open, since our Hamiltonian, Eq. (T4]), does not include 
the effect of the scattering potentials due to the holes on the external shell on the nanotube 
electrons. In principle, the scattering of the electrons by these holes should be included. Our 
approximation, which neglects these effects, is based on two equally important observations. 

The first observation is that there is experimental evidence, to suggest that electron 
transport through a peapod in a two-terminal conduction geometry, and in presence of a 
gate electrode, is ballistic over the entire sample (whose length is around ~ 400 nm). In 
other words, the holes on the nanotube shell do not hinder ballistic electron transport. 

The second observation is the following: the inclusion in the model of any topological 
defects (holes) in the tube shell would certainly affect the spectrum of the electrons in the 
nanotube and the level spacing. However this would only result in a renormalization of the 
coupling constant A, which in our model is a fitting parameter. The anomalous temperature 
dependence of conductance induced by vibrational dynamics of encapsulated fullerenes will 
not be affected. 
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Finally, we have also showed that the description of the hybridization-induced scattering 
potential as a long-range potential is not adequate to describe polaronic effect. 
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